Nonequilibrium transport through a singlet-triplet Anderson impurity 
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We study the Anderson model in which a configuration with a doublet is hybridized with another 
with a singlet and a triplet. We calculate the conductance through the system as a function of 
temperature and bias voltage, near the quantum critical line for which the system is exactly solvable. 
The results explain recent transport measurements in a single-molecule quantum dot. 
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INTRODUCTION 

When the wave function is forced to evolve continu- 
ously between two competing ground states, a quantum 
phase transition (QPT) takes place between these two 
states At the transition, the length scale of quantum fluc- 
tuations becomes infinite, and exotic states of condensed 
matter are expected Recently, several remarkable 
features of a (QPT) have been observed in Cgo quantum 
dots (QD's) with even occupancy inserted in a nanoscale 
constriction ■ In the last years there has been a great 
interest in systems of quantum dots (QD's) because of 
the possible technological applications, and also because 
they constitute ideal systems with a single magnetic im- 
purity in which several parameters can be tuned. 

When the QD has an odd number of electrons and the 
Coulomb repulsion U is large enough, the conductance at 
zero bias is increased below a characteristic Kondo tem- 
perature Tk as a consequence of the Kondo effect. This 
is a usual feature of single-electron transistors built with 
semiconductor QD's or single molecules 0,01 and is well 
understood in terms of the simplest ordinary Anderson 
model (0AM) in which a configuration with a doublet is 
hybridized with a singlet. In a dot with an even num- 
ber of electrons, there are two competing states for the 
ground state configuration: a singlet in which two par- 
ticles occupy the lowest level (say |00) — s|sj|0)) and 
a triplet in which one electron is promoted to the next 
level and coupled ferromagnetically (|11) = p|s||0) and 
its SU(2) partners) due to the strong Hund couphng Q. 
The simplest Anderson model which describes the system 
mixes these four states with a doublet (|a) = sj^jO)) by 
promoting a particle (electron or hole) to one of the leads 
i,0|. This is the singlet-triplet Anderson model (STAM) 
which had been used to describe valence fluctuating Tm 
impurities in a cubic environment 0. 

This model has a quantum phase transition (QPT) 
from a singlet to a doublet ground state as the energy 
of the triplet is decreased [a, When the triplet is 
well below the other states, there is a partial screening 



of the spin 1 that explains the zero bias Kondo peak ob- 
served experimentally in this situation On the 
other (singlet )side of the transition there is a dip in the 
conductance ^, [3| that has also been observed experi- 
mentally (1, [7|. 

In this paper, we show that the experimental obser- 
vations in Geo QD's can be understood in terms of this 
model and its QPT. The differential conductance dl/dV 
as a function of temperature T and bias voltage V has 
been measured at both sides of the transition 0]. On 
the singlet side of the transition, a dip in the conduc- 
tance at y = is observed in agreement with theoretical 
expectations 0] as well as non-equilibrium measure- 
ments performed in carbon nanotubes 0]. On the other 
side of the transition, dl/dV as a function of V shows 
a structure with three peaks that has not been quan- 
titatively explained yet. As the temperature T is de- 
creased, the zero bias conductance G{T) first increases, 
then shows a shoulder or a plateau and then increases 
again. As stressed by the authors, this behavior is still 
not understood [5] . The last two figures of our paper are 
the theoretical counterpart of these experimental results. 
The main qualitatively features are reproduced. A more 
quantitative agreement would require a fine tuning of the 
parameters which is beyond the scope of this work. We 
provide an interpretation for the observed behavior. 



THE MODEL 

The STAM assumes infinite U and contains two neigh- 
boring charge configurations in the QD: c?" and c?"+^ with 
n and n + 1 particles. Without loss of generality we can 
assume n to be even, performing an electron-hole trans- 
formation if necessary. Then, the configuration con- 
tains a singlet \SM) = |00) , where S is the spin and M its 
projection, and a triplet |1M), (M = —1, or 1), while 
the d""*"^ configuration consists of a doublet denoted by 
its spin 1/2 projection \a). We introduce the following 
creation operators for a particle in the dot 11 1 
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k)(oo|, 



< = -(|T)(10| + ^/2U)(l-l|)/^/3, 
(U>(10| + V2| T>(ll|)/V3. 



4i 



(1) 



The operators dl^ and hybridize via matrix elements 
VJi. and Vjj, with the conduction states c^fco- of two con- 
ducting leads V = L (left) or R (right) that transport the 
current through the QD, leading to the Hamiltonian 



H = Esmm + EtY.\lM){lM\+EdY,\'^){^\ 

M a 

+ E \y^kdL + VlAaY^u. + H.c 



(2) 



We assume ^^^.V^j, — ^Lk^Rk^ that only one con- 



duction channel ~ X^i, 



(ry = s or t) hybridizes 



with the dot states. In general, also the orthogonal linear 
combination of c^ka- plays a role and "screens" the re- 
maining doublet ground state when the localized triplet 
is well below the singlet, leading to a singlet ground state 
13, l^- However, the energy scale involved in this second 
screening T* (which depends exponentially on a small 
coupling constant [isj ') might be very small. As discussed 
in Ref. [iSj, this is likely the case of previous experiments 
0) as well as those in Ceo QD's the theory in 
the general case [H, [l^l predicts that the zero-bias con- 
ductance G(T) should decrease at very low temperatures 
and dl/dV should also decrease for the smallest applied 
bias voltages V in contrast to the observations. This in- 
dicates that T* is smaller than the smallest temperature 
in the experiments. 

The STAM with only one conduction channel, also de- 
scribes the mixing between the low lying states of the 
4f^^ and 4i^^ configurations in a cubic crystal field [1]. 
For Et +00, the model reduces to the 0AM. For 
Es — > -|-c», the model describes valence fluctuations be- 
tween two magnetic configurations 14]. In both limits, 
for constant density of conduction states and hybridiza- 
tions, the model is exactly solvable (by the Bethe ansatz) 
and the ground state is a singlet (doublet) in the first 
(second) case Thus, the model has a QPT as a func- 
tion of Eg — Et- The position of the transition depends on 
the other parameters of the model, leading to a quantum 
critical surface that can be determined calculating the 
magnetic susceptibility at T ^ using numerical renor- 
malization group (NRG) Q. However, li \Vt\^ = ■i\Vs\^ , 
the transition takes place exactly at Eg — Et ~ inde- 
pendently of the value of Ed- In addition, along this line, 
the model can be mapped into an 0AM plus a free spin 
1/2 % We will use these results to control the distance 
to the QCP. 



APPROXIMATIONS AND EQUATION FOR THE 
CURRENT 



As discussed for example in Ref. Il5|, the calculation of 
non-equilibrium properties of a strongly correlated sys- 
tem is a particular challenge for theory. A recent ex- 
tension of the numerical renormalization group (NRG) 
for the non-equilibrium case seems promising, but is not 
fully developed yet [3]. For the 0AM, the non-crossing 
approximation (NCA) llTII and renormalized perturba- 
tion theory in U (isl Il9j have been useful. However, 
the latter is very difficult to extend to the STAM. The 
so called poor man's scaling works well when either eV 
or the magnetic field energy is larger than kTx [lOl and 
has been successfully extended for a model similar to the 
STAM on the singlet side of the QPT 0. However, this 
method ceases to be valid near the quantum critical sur- 
face for small V . Instead, the NCA can be extended to 
more general Anderson models, like tose appropriate for 
Ce compounds 21|, Co impurities on Ag and Cu 12211 or 
systems of two quantum dots out of equilibrium [23[ . 

In this work, we extend the NCA to the STAM out 
of equilibrium. We introduce auxiliary bosons, one 
for the singlet state and three for the triplets, and 
auxiliary fermions for the doublet, in analogy to the 
SU(N)xSU(M) generalization of the Anderson model 
[2^ 1 . The spectral densities of the operators d\^^ defined 
by Eqs. H]) for given spin, p^(w) and p^(w), are deter- 
mined by convolutions from those of the auxiliary parti- 
cles. The current is given by [l^l 



1 = ^ /d^[rv^(c.) + rv*(c.)][/i(c.)-/«(c.)], (3) 



where = T\ + Tl with = \VI\\^5{uj - ek) as- 

sumed independent of u within a bandwidth D and zero 
elsewhere, the asymmetry parameter A = 4r^r^/(r^ -t- 
r^)^ (independent of 77), and /j/(w) is the Fermi function 
with the chemical potential of the corresponding lead. 



RESULTS 

For the numerical evaluations, we consider the case in 
which the ground state configuration has an even number 
of particles {Es,Et < Ed). We choose Ed = + Mfl) = 
(without loss of generality) and take a band width 
D = lOr, l/2r, and F* = 3/2F, where F will 

be our unit of energy. The choice F*/F* — 3 allows 
us to have an accurate control of the distance to the 
QPT, while the main features of the results depend on 
this distance and not on the specific choice of parame- 
ters. For F7F" = 3, if in addition Eg = Et, the STAM 
can be mapped into an 0AM plus a free spin 1/2 (Ref. 
;8]). The 0AM has total coupling F"^" = F and in- 
verted charge transfer energy and chemical potentials: 
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FIG. 1: (Color online) singlet (left) and triplet (right) contri- 
butions to the dot spectral density for Et — —3, Es — —2.9 
and different temperatures. The flatter curve corresponds to 
the highest temperature. F = 1 is the unit of energy. 
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FIG. 2: Differential conductance as a function of bias voltage 
for several temperatures and (a) Et = —2,Es = —2.03, (b) 
Et = -3,Es = -3.1. 



^OAM _ ^OAM ^ Es~ Ed, /i?^^ = -M^- Using this 
mapping it can be shown that the densities of both mod- 
els are related by p^i^j) = p^iuj) = /9°^^(-tj)/2 and 
the absolute value of the current is the same. From the 
structure of the corresponding NCA equations for both 
models, we realize that the NCA satisfies these equalities. 
This has also been verified numerically. Furthermore, for 
any Ed the QPT takes place at E^ = Et when r'/r" = 3 
Q . This defines a quantum critical line and then Es — Et 
controls the distance to this line. 

As explained above, for F = (/xl = /ir = 0), when 
Eg — Et, the singlet and triplet parts of the spectral 
density of the dot 2p^{Lu) and 2p^^{u}) coincide with the 
mirror image of the localized spectral density alre ady re- 
ported for the Anderson model (Fig. 5 of Ref. 17[ for 
Et = — 2r) and exhibit the usual Kondo resonance at 
the Fermi level. The half width of this peak allows to 
define a Kondo temperature Tk- How do the spectral 
densities evolve as one moves away from the quantum 
critical surface? Decreasing Eg (on the singlet side of 
the QPT) Pdi^) displaces to positive frequencies, while 
/9^(a;) decreases and displaces its weight to negative fre- 
quencies. As a consequence, a gap opens in the sum 
r''p^(u;) -fr*/9^(a;) entering Eq. ([3]) at low temperatures. 
This situation has already been studied previously and 
interpreted as the low temperature part of a two-stage 
Kondo effect [1]. 

The spectral densities for Es > Et are shown in Fig. [T] 
In contrast to the previous case, p^(tf) remains peaked 
at the Fermi energy at low temperatures. This is a con- 
sequence of the partial Kondo effect, by which the spin 
1 at the dot forms a ground state doublet with the con- 
duction electrons of both leads [IJ]. The singlet part 



of the density p^(w) displaces to negative frequencies in 
this case. Therefore a pseudogap appears in the sum 
Vp^^uj) + T*pIi{uj), but at finite frequencies, in contrast 
to the gap at the Fermi level that develops at the singlet 
side of the QPT. Note that at high temperatures both 
densities are quite similar and the differentiation between 
Pdi^) and Pdi^) develops at a characteristic temperature 
of the order of a fraction of Es — Et- 

The equilibrium (V = 0) conductance as a function of 
temperature G{T) on the singlet side of the QPT, near 
the quantum critical line (not shown) shows a maximum 
at a finite temperature and agrees with previous results 
using NRG Q and with experiment Q- In particular, the 
increase and decrease of G(T) from its maximum value 
are logarithmic to a good degree of accuracy. 

The differential conductance dl/dV as a function of 
bias voltage on the singlet side of the transition is dis- 
played in Fig. [2] (a) . We have applied the voltage sym- 
metrically (hl = —PR — As the temperature is 
lowered, a dip develops at small voltages. The half width 
of the dip is of the order of Et — Es- These results are in 
good agreement with the experimental ones (Fig. 4 c of 
Ref. [2|). For other parameters, in particular larger val- 
ues of Et — Es and Ed — Es (less valence fluctuations) we 
obtain curves that look similar to those reported in car- 
bon nanotubes, with a flat bottom at low temperatures, 
explained using poor man's scaling An example is 
illustrated in Fig. [2] (b) 

The most novel experimental results are those on the 
"triplet side" of the transition {Es > Et). As the tem- 
perature is decreased, G(T) increases until it reaches a 
plateau at the characteristic energy Es — Et and then it 
continues to increase j2| . This is in general, the behavior 
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FIG. 3: (Color online) Zero bias conductance G{T) as a func- 
tion of temperature (full line) and contributions from the 
singlet (dashed dot dot line) and triplet (dashed line) for 
Et = -2>,Es = "2.9. The inset shows G{T) for Et = -3 
and several values oi 5 = Es — Et (increasing from from top 
to bottom). The circles correspond to Eq. Q. Go = 2e^ A/h. 



that we obtain, as displayed in Fig. [S] As the partial 
contribution to the conductance for each density reveals, 
the plateau is due to the contribution of the singlet part 
of the dot spectral density, namely p^{oj) which is peaked 
at —{Es — Et) (see Fig. [T|). This plateau has not been 
observed in previous calculations using NRG fi]. We be- 
lieve that the reason for this is the lack of resolution of 
NRG to describe peaks in the spectral density out of the 
Fermi energy, like that of p^{(^)- A clear manifestation 
of this is a system in which the Kondo peak is split in 
two out of the Fermi level [26|. In fact, while NRG uses a 
logarithmic frequency mesh centered at the Fermi energy, 
we find that in order to obtain enough accuracy in the 
convolutions that define p^(ti') and p^(cli), it is necessary 
to use two different dense logarithmic meshes, centered 
at the corresponding peaks. 

The inset of Fig. [3] shows the evolution of the equi- 
librium conductance as the system is displaced from the 
quantum critical line to the triplet region. At this line, 
G(T) is the same as the corresponding result for the 0AM 
obtained using the mapping mentioned above. In this 
case, a fit of G(T) using the empirical curve derived by 
fitting results of the NRG for a spin 1/2 
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with s — 0.22 works very well. We could not find other 
regions in which similar empirical curves would fit well a 
large portion of the curve. When Eg increases, the degen- 
eracy in the ground state is reduced and G(T) decreases 



0.1 



; 0.08 - 



> 

■a 

^ 0.06 



0.04 





1 


1 


' 1 ' 

- 1=0.01 r 








-- 1=0.005 r 








T— n nnn^ r 

1 — U.UUUJ i 


- 

- 






- 












'■'/ 






/ — \ ■■ '/ 


\\ ■. / \ 




/> • ■ ^ y / 

1 , 1 


\ ^/ ^ ^ 

1 , 1 



-0.2 



-0.1 





eV/r 



0.1 



0.2 



FIG. 4: Differential conductance as a function of bias voltage 
for Et — —3, Es — —2.9 and several temperatures. 



in the whole range of temperatures. However, in contrast 
to the case Eg < Et, the conductance at zero tempera- 
ture retains values near to the ideal one Go = 2e^A/h 
due to the partial Kondo effect. 

In Fig. [4] we show the non-equilibrium differential con- 
ductance dl/dV on the triplet side at several tempera- 
tures. At low temperatures, there is one peak centered at 
V = and other two centered at V = ±Vm- This three- 
peak structure also agrees qualitatively well with experi- 
ment 0. Actually in the later, the peak at = —Vm is 
higher than that at = Vm, while our results are even 
in V due to our assumption of a symmetric voltage drop. 
Changing this we can easily control the relative height of 
both peaks. In any case, the main point is the existence 
of these three peaks, which is also a consequence of the 
particular structure of the spectral densities and 
p^j^iuj) at equilibrium and low temperatures (see Fig. [T]). 
The peak at F = is due to the pli{uj) which is peaked at 
the Fermi energy. When e\V\ reaches energies at which 
p^{oj) is peaked, this density starts to contribute signifi- 
cantly to the current / [see Eq. ([3|)] and dl/dV increases. 
Finally, for large e\V\, I tends to saturate and dl/dV de- 
creases again. 

The effect of temperature is to broaden the peaks and 
at high enough temperatures only one broad peak in 
dl/dV is present, in agreement with experiment. For dif- 
ferent parameters as those of Fig. 1, this broadening of 
the spectral densities is already important at = T = 
and increases with |V^| in such a way that even at T = 
only one peak in dl/dV is present. This happens for ex- 
ample at Et = -2r, Es- Et = O.Oar, for which dl/dV 
(not shown) has a monotonic behavior for positive V , 
displaying only a shoulder for eV ~ ±(£^5 — Et). 

In summary, the transport properties recently ob- 
served near the singlet-triplet quantum phase transition 
in quantum dots can be explained in the framework of the 
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singlet-triplet Anderson model with one channel per con- 
duction lead, out of equilibrium, using the non-crossing 
approximation. We made use of exact results to control 
the distance to the quantum critical point. The differen- 
tial conductance dl/dV as a function of bias voltage is 
markedly different at both sides of the transition show- 
ing a dip (peak) at small voltages on the singlet (triplet) 
side and often a three-peak structure on the triplet side. 
The zero bias conductance G{T) as a function of temper- 
ature displays a plateau due to the contribution of the 
excitations from the localized singlet to the doublet. 
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